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Abstract 

We present a fairly new and comprehensive approach to the study of stationary flows of the 
Korteweg-de Vries hierarchy. They are obtained by means of a double restriction process from 
a dynamical system in an infinite number of variables. This process naturally provides us with 
a Lax representation of the flows, which is used to find their bi-Hamiltonian formulation. Then 
we prove the separability of these flows making use of their bi-Hamiltonian structure, and we 
show that the variables of separation are supplied by the Poisson pair. 



1 Introduction 



In this paper we present a comprehensive approach to the stationary flows of the Korteweg-de Vries 
(KdV) equation. This is a quite classical subject, and lies at the very heart of the modern theory 
of integrable systems in finite and infinite dimensions. Its formulation can be traced back at least 
to H 22]. Then, together with its generalizations, it has been subject of intensive studies. See, for 
example, |7|, 12, |l3[ |b|, [2^, 28] and references therein. It was realized that such systems were among 



the prototypes of the class of Algebraically Completely Integrable Hamiltonian Systems |2], p9[ , and, 
in particular, it has been proven that the classical Jacobi formulas for projective embeddings of 
hyperelliptic Jacobian varieties found a very natural realization in such problems. 



Hamiltonian aspects of those flows were studied in a number of papers. To cite a few || 14, 23]. 
Two main approaches emerged. In the first one, a Hamiltonian structure was given to the stationary 
KdV flows, looking at their variational properties ||. Namely, such flows were regarded as classical 
Euler-Lagrange equations associated to suitable Lagrangian and Hamiltonian densities. In the 
second one (see, e.g., ]l3|), a set of suitable canonical coordinates was introduced on the stationary 
manifolds, somehow dictated by the algebro-geometrical structure of the problem, and the fact 
that the flows were indeed Hamiltonian was verified at a later stage, via direct computation. In 
the same way, a set of action-angle variables was found. 

After the discovery of the bi-Hamiltonian structure of the 1 + 1 dimensional KdV hierarchy, 
the obvious problem of finding a corresponding bi-Hamiltonian structure for its stationary flows 



was studied. A solution for this was found 0, p2l 37] with an ingenious study mainly based on the 



recursion relations found by Alber p], the introduction of special sets of coordinates, and the use 
of the Miura transformation. 

The aim of this paper is to give a rather new and systematic perspective to this circle of ideas, 
focusing our attention on the bi-Hamiltonian aspects of the problem. The next section is devoted to 
a quick description, in a simple example, of the stationary reductions of KdV and to an illustration 
of the main properties of these systems. Then, we will present the plan and the important points 
of the paper. 



2 A Preliminary View 

In this section we recall some known facts about the KdV hierarchy and its stationary reductions, 
and we present the problems to be tackled in the next sections. 

The KdV equation is the most famous example in the class of the so-called integrable nonlinear 
PDEs. It possesses a number of remarkable properties, in particular: 

1. It has an infinite sequence of integrals of motion; 

2. It admits a Lax representation; 

3. It is a bi-Hamiltonian system; 

4. The integrals of motion are the coefficients of a Casimir of the Poisson pencil, so that the 
KdV equation can be seen as a Gel'fand-Zakharevich system (see below). 

There are of course relations between these properties. For example, the conserved densities are 
the residues of the fractional powers of the Lax operator. Also, they can be extracted from the 
bi-Hamiltonian structure and shown to commute with respect to both Poisson brackets. The 
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associated vector fields form the KdV hierarchy, whose first members are 
du du 1 

^r = ^' -ET = li u xxx - 6uu x ) (KdV equation) 

OU 1 2 

wr = —{u xxxxx - 10uu xxx - 20u x u xx + 30-u u x ) . 
OT5 lb 

The KdV hierarchy can be used to find finite-dimensional reductions for the KdV equation, giving 
rise to explicit solutions. Indeed, the set of singular points of a (fixed) vector field of the hierarchy 
is a finite-dimensional manifold which is invariant under the flows of the other vector fields, due 
to the commutativity property. The (finite-dimensional) systems obtained by restricting the KdV 
hierarchy to such invariant manifolds are called the stationary reductions of KdV . 

Let us consider explicitly the reduction KdVs corresponding to the third vector field of the 
hierarchy. The set of zeroes is given by 

^xxxxx ^-Q^^xxx ^^^x^xx 30zi ^J>x — 5 ^2.2^ 

and its dimension is 5, since we can use the value of u, u x , u xx , u xxx , and u xxxx at a fixed point xq 
(i.e., the Cauchy data) as global coordinates. For the sake of simplicity we put 

u = u(x ), ui = u x (x ), u 2 = u xx (x ), u 3 = u xxx (x ), u A = u xxxx (x ) . (2.3) 



In order to compute the reduced equations of the first flow of Q2.1| ), we have to derive it with respect 
to x, and to use the constraint fl2.2|) and its differential consequences to eliminate all the derivatives 
of order higher than 4. We obtain the equations 



duo du\ du2 du 3 

Ot\ Oti Oti Ot\ 

dui 2 

— — = WuqU 3 + 20liiU2 — 30uo Ml . 
Oti 



(2.4) 



In the same way, for the KdV equation we get 

du 1, a s 
— = -(u 3 -6u 0Ul ) 

^ = i(u 4 - 6u u 2 - 6ui 2 ) 
du 1 

— — = -{^uqu 3 + 2uiu 2 - 30u 2- "i) (2.5) 
dt 3 4 

— ^ = -{Auqua + 6u±u 3 + 2u2 2 — 30uo 2 U2 — GOuoui 2 ) 
oh 4 

dui 1 2 1 3 

— — = -(10ui«4 + lOtio u% + 10n 2 U3 — 100uqUiU2 — 60u\ — 120uo «i) 
dt 3 4 

As far as the restrictions of the other flows are concerned, it can be shown that they are linear 



combination of (2.4) and (|2.5| ). 

It is not surprising that the above mentioned properties of the KdV equation hold also for the 
KdVs system. However, to the best of our knowledge, it has not been made completely clear the 
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way which these properties pass from the KdV hierarchy to its stationary reductions, especially as 
far as the bi-Hamiltonian structure is concerned. In any case, one can check that the functions 



Hq = ( — 2li2«4 + 6uo 2 U4 + U3 2 — 12uqUiU3 + 16uqU2 2 + \2u\ 2 U2 — 60uo 3 U2 + 36uo 5 ) 

16 



Hi 



-{2uQUi — 2u\uz + U2 — 20uo 112 + 15uq 



(2.6) 



H2 = U4 — WuqU2 — 5u\ 2 + IOuq" 



are integrals of motion for (2.4) and ( |2.5| ). These systems have also a Lax formulation, i.e., they 
can be written as 



8L 

dti 



1,3, 



(2.7) 



where the Lax matrix L depends on a parameter A, and is given by 



L 



16 



4«iA + 113 — Ququ\ 



16A 3 + 8n A 2 + 2A(u 2 - u 2 )+ 
V 1i4 — SUQU2 — 6ui 2 + 6uq 3 



16A 2 - 8u A + 6u 2 - 2u 2 ^ 
— 4ttiA — 113 + GuqUi 

J 



(2.1 



The matrices Ai can be easily contructed from L (see Section ||). 

Finally, there are two compatible Poisson structures giving a (bi)-Hamiltonian formulation of 
the KdVs systems. The corresponding Poisson tensors are 



and 
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If we call X\ and X3 the vector fields of KdVs, then the following relations hold: 



P dH 2 = 
X l = P dH t = P 1 dH 2 
X 3 = P dH = P 1 dH 1 

PidH = . 

They can be collected in the statement that the function H(\) := H2X 2 + H\ \ + Hq is a Casimir 
of the Poisson pencil P\ := Pi — APo 5 that is, 

P x dH(X) = . (2.10) 

Therefore, X\ and X% are the bi-Hamiltonian vector fields associated with a polynomial Casimir 



of a Poisson pencil of maximal rank. In a word, they are Gel'fand-Zakharevich (GZ) systems [19]. 

The importance of the stationary reductions of the KdV hierarchy (and, more generally, of the 
stationary reductions of the Gel'fand-Dickey hierarchies) lies in the fact that the reduced equations 
can be solved by means of the classical method of separation of variables. This was noticed in the 
early works on the subject. It is also known how to construct the variables of separation starting 
from the Lax matrix. We will show that the separability of these systems is a particular instance 
of a general result, which is valid for quite a wide class of bi-Hamiltonian manifolds. 

In Section ^ we give a rather unconventional presentation of the stationary reductions of KdV. 
Our priviledged starting point is a picture of the KP hierarchy as a system of ordinary differential 
equations, called the Central System (CS) in j|, |l5f . Starting from there, by means of a double 
reduction process we can describe quite explicitly the stationary reductions of KdV, and, in Section 
H], we are able to give a Lax representation of these systems, with a Lax matrix depending poly- 
nomially on a parameter. This representation is used to show that the flows are bi-Hamiltonian. 
This is done in two steps. First, in Section ^ we recall the bi-Hamiltonian structure on matrix poly- 
nomials and show that the Hamiltonian vector fields (with respect to the Poisson pencil) admit 
a Lax formulation. This property is conserved after a suitable bi-Hamiltonian reduction process. 
Then, in Section ||, we identify the phase space of the stationary reductions of KdV with a re- 
duced bi-Hamiltonian manifold, and we show that they are GZ systems. In Section ^ we state 
(referring to [16] for a more detailed discussion) a theorem ensuring that, under some additional 



assumptions, the GZ systems are separable in coordinates that are naturally associated with the 
bi-Hamiltonian structure. Finally, in Section |8] we show that this theorem can be applied to the 
stationary reductions of KdV, and that the variables of separation can be constructed algebraically. 

Summing up, in this paper we present a somewhat self-contained approach to the study of 
stationary flows of KdV, and we use them as a laboratory to test ideas of the bi-Hamiltonian 
geometry, from the GZ theory to the separation of variables. In our opinion, such a set up provides 
a comprehensive formulation of results which, although for the most part already available in the 
literature, would perhaps acquire a deeper meaning under this perspective. 



3 KdV Stationary Reductions 

In this section we give a self-contained presentation of the stationary reductions of the KdV hierar- 
chy, using the formalism developed in || |15|] for the KP theory. Our starting point is the Central 
System (CS), a family of dynamical systems with N x N degrees of freedom. A first (stationary) 
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reduction gives rise to the CS2 hierarchy, with N degrees of freedom. Then a further restriction 
leads to finite-dimensional systems that coincide with the stationary reductions of KdV. 

We consider the space TL of sequences {-H^}fc>i of Laurent series having the form = 
z k + Yli>i Hfz~\ where Hf are (complex) scalars that play the role of coordinates on TL. On such 
phase space TL we define a family of vector fields as follows. We associate with a point {H^}k>i 
in TL the linear span H+ = {H^\H^,H^\...), where # (0) = 1. The defining equation for the 
j-th vector field Xj of the family, to be referred to as the Central System (CS), is the invariance 
relation 

J- + «»)« + C« + . (3.1) 
This relation is equivalent to the (explicit) equations 

= _ H U) H (t) + H U+k) + J2 fl-f fl"(*-0 + H?H^- l \ k > 1. (3.2) 
3 1=1 1=1 



Remark 3.1 From (|3.2|) it is evident that the exactness property 



dt k H ~ dtj [6 - 6) 



holds. Moreover, it can be shown [H] that the flows of the CS commute. 



Remark 3.2 There is a very tight relation between the CS and the linear flows on the Sato 



Grassmannian [34, 35, 11]. This relation is discussed in [15], where the classical result of Sato on 
the linearization of the KP hierarchy is recovered from the point of view of the bi-Hamiltonian 
geometry. 

Since the CS is a family of commuting vector fields, we can reduce it in many different ways. By 
means of a suitable combination of such reduction processes, the so-called fractional KdV hierarchies 
were obtained in || . Now we will show how the stationary reductions of KdV can be derived from 
the CS. The commutativity of the flows implies that the set Z2 of zeroes of the vector field X2, 
defined by the quadratic equations 

H (k+2) _ H (k) H (2) + HfH^ + ffftfW + Hi = 0, (3.4) 
l=i 

is an invariant submanifold for CS. Moreover, on Z2 we have 

BHi2) dH ' f ' 0, (3.5) 



dtj dt2 

due to the exactness property fl3.3|) . Therefore, the manifold Z2 is foliated by invariant submanifolds 
defined by the equation = constant. Among all these leaves, the submanifold 52 defined by 
the simple constraint 

= z 2 (3.6) 
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is particularly relevant. At the points of £2 equation (|3.4| ) takes the form 

H (k+2) = z 2 H {k) _ H k R (l) _ H k ? (37) 

and allows us to recursively compute the Laurent coefficients of H^ k \ for k > 2, in terms of the 
coefficients of h := H^K Hence, £2 is parametrized by the coefficients {/iz};>i of h. Equation ( |3.7| ) 
also shows that z 2 (H + ) C H+, so that on £2 the elements {z 2 i , z 2 ^h}j>Q form a basis in H + . Thus, 
we have that 

H {k) = Pk (z 2 ) + q k (z 2 )h onS 2 , (3.8) 

where p k and q k are polynomials. This can also be seen directly from equation (3.7). Moreover, 
there is only one Laurent series of the previous form satisfying the asymptotic condition = 
z k + 0(z^ 1 ) as z — > 00. 

Definition 3.3 (see [§]). The restriction of CS to the invariant submanifold S 2 is called the CS 2 
hierarchy. 

The restricted vector fields are given by 

™ = _ hH (i) + #Cj+D + £ hlH U-D + Hi j > h (3 .9) 
■? i=i 

where the must be written in terms of h according to (p.Sj). If we denote with H- the span 
of the negative powers of z, and with 7r_ the projection on H- according to the decomposition 
H+ H— , then the equations ( |3.9| ) can be written in the more compact form 

^ = -7T_ { qj h 2 ) . (3.10) 

Notice that _ff( 2fc ) = z 2k , so that ( |3.3| ) implies that the even flows of CS2 are trivial. 

The finite-dimensional systems that are the main subject of this paper are those obtained by 
restricting the CS2 flows to the manifold of zeroes of the {2g + l)-st vector field of CS2. We will 
call such systems the KdV2 ff +i systems, since they are (equivalent to) the stationary reductions of 
the KdV hierarchy, as we are going to show at the end of this section. The constraint which defines 
the phase space M. 2g +i of the KdV2 9 +i system is 

r>h 

«P- = (Q2 9+ ih 2 ) = . (3.11) 

Ot 2g+ l 

A direct inspection shows that this constraint gives all the coefficients of h in terms of the first 
(2g+l), i.e., hi, ... , h 2g+ %. In other words, the dimension of the phase space of the KdV2 9 +i system 
equals 2g + 1. The equations are given by the first 2g + 1 components of ( 3.10| ), after substituting 
the constraints ( |3.11| ). In the case of KdVs there are two independent vector fields: 

dh\ , dhi 

—± = -2h 2 ^ = -2/14 + 2hi h 2 

oti at 3 

dho , , 9 dho , , 9 , o 

2 = _2h 3 -hi 2 — ± = -2h 5 + h 2 2 + V 

oti at 3 

^ = -2hih 2 - 2/14 ^ = -2/11/14 + 4/ii 2 /i 2 -2h 3 h 2 (3.12) 

oti at 3 

^ = -2/15 - h 2 2 - 2h 1 h 3 ^ = -2/13 2 - 2h 2 h 4 + 2hih 2 2 + hi A + hi 2 h 3 

oti at 3 

Bh Bh 

— - = -Ah 3 h 2 + 2hi 2 h 2 - 4/11/14 = 2/ii 2 /i 4 - 4/13/14 + 2/ii 3 /i 2 

ail ot 3 
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These are, up to the coordinate change (|3.14|) , the equations (|2.4|) and (|2.5|) . 
We remark that along the flows of KdV2 9 +i the relations fl3.3|) take the form 



9ff(29+l) 

-gjj- = 0. (3.13) 

showing that all the coefficients of H^ 2g+1 ^ are integrals of motion. Therefore our presentation of 
the KdV stationary reductions carries directly the conserved quantities of the flows. Moreover, in 
the next section we will show that the Lax representation also arises in a natural way. We end this 
section with the following: 

Remark 3.4 The usual KdV hierarchy in 1 + 1 dimensions is described in j|] as a projection of 
CS2 along the integral curves of the first vector field of the hierarchy, 

_ = -h 2 + z 2 + 2h x . 

Oti 

Indeed, if we put x = t± and u = 2h±, then the previous equation takes the form h x + h 2 = u + z 2 
and allows us to write the hj as polynomials in u and its x-derivatives: 



hi = 


\u 




h 2 = 


~ ~g U X 




h 3 = 


g (jJ'XX U 


) 


/14 = 


16 i^xxx ~ 


- 4:UU X 


h 5 = 


^(tlxxxx 


6uu xx 



(3.14) 

hu 2 x + 2u 3 ) 



Thus equations ( |3.9| ) become partial differential equations for the variable u, and are the KdV 
hierarchy. But we can also use the system (3.14) to recover CS2 from the KdV hierarchy, so that 



we can pass back and forth from one hierarchy to the other. This shows that the KdV2 9 +i systems 
that we have introduced coincide with the usual stationary reductions of KdV. The first (2g + 1) 
equations of the system ( 3.14 ) represent the change between our coordinates (hi, . . . , fog+i) and 



the ones usually considered in the literature, namely (u, u x , . . . , u^ 29+l ^). 



4 The Lax Representation 

In this section we show that there is a quite natural (Zakharov-Shabat) zero-curvature representa- 
tion for the CS2 system, entailing a Lax representation for the KdV2 5 +i hierarchy. 

We know from the previous section that in the CS2 theory every element in H + can be written 
as a linear combination of 1 and h with coefficients that are polynomials in A := 
point of the manifold 52 (that is, to each series h = 



. Then to each 
^ ;>1 hiz~ L ) we can associate a family of 



2x2 matrices V^^(A) depending polynomially on A, defined by the relation 



1 




1 


h 




h 



(4.1) 
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Since the even flows are trivial, we will be interested only in the matrices of odd index. The first 
three of them are given by 



v( 3 ) 



v( 5 ) 



1 

\ + 2h 1 

— h 2 X — Ha + h\h2 
X 3 + h 1 \ 2 + h 3 X+ 



—hi A — hi 

A 2 + h\\ + 2/i 3 - hi 2 h 2 

A 2 - h x X - h 3 + h x 2 
/12A — h\h2 + hi 



2h 5 - 2/ii /i 3 - h 2 2 + hi 3 
The commutativity of the flows and the "abelian" zero-curvature relation 



imply that 



_d_ 



_d_ 



yU) + v (fc) ,V (i) 



1 



0. 



(4.2) 



Since the entries of the matrix appearing in the previous equation are polynomials in A, and the 
elements {X J , X J h}j>o are linearly independent in H+, it follows that the zero-curvature relations 



Of 



dt k 







(4.3) 



hold. 

If we restrict to the set .M23+1 of the stationary points of the (2g + l)-st vector field of CS2, 
the zero-curvature representation naturally gives rise to Lax equations for the matrix \A 2ff+1 ) , 



_d_ 

dtl 



y{2g+l) = r v (fc) )V (2g+l) 



(4.4) 



The following proposition will be useful in Section ^, and shows that these Lax equations faithfully 
represent the KdV2 3 +i system. 



Proposition 4.1 The matrices \A 2fc+1 ) of the CS2 hierarchy have the following properties: 



1. The matrix \/( 2k+1 ^ depends only on (hi, . . . , /i2fc+i)j <ind the map (hi, 
is infective; 

2. The trace of V( 2fe+1 ) is zero: 

3. For i < k one has 



V (2fe+1) 



V (2i+1) = (A i-fe v (2fc+l) )+ _ a 



ik 





1 



(4.5) 



where (•)+ denotes the projection on the nonnegative powers of X, and is the entry (1,2) 



of the coefficient of X 



k-i-l 



III 



V (2fe+1)_ 



Proof. First of all we observe that, almost by definition, 
y&k+l) = 



P2k+1 



Q2k+1 



, ' , n 2k+l 1*1=1 h 2iq2k-2l+l 

1^1=1 h2lP2k-2l+l + til 



(4.6) 
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Then, we notice that equation (|3.7f ) implies the recursion formulas 



(4.7) 



Thus, by induction, we obtain 



Ptt+i = - E ^~ lAfc ~' > **n = A " - E a?- 1 **-' • 

J=l 1=1 

In order to express the coefficients H^ 1-1 and Hf 1 ^ 1 in terms of the fy, we use the identity 

k+l-3 

Tj2k-1 _ rrl \ rrl rr2fc-2i-l 

n l — n 2k+l-2 ~ 2-^1 n l+2i~2 n l > 
i=l 

which can be proved by induction on using again (|3.7|). In particular, we have 



(4- 



fc-2 



H 



2fc-l 



#2fc-l - E H li-t H l 



2k-2i-\ 



i=l 
fe-1 



if. 



2fc-l 



#2fc ~ E ^2i-^l 



2k-2i-l 



(4.9) 
(4.10) 



i=l 



This allows us to control the appearance of the h{ = Hf in P21+1 and q2i+i, and leads to the proof 
of the first assertion. 

The second statement tantamounts to P2fc+i + X)f=i ^2Z<?2fc-2Z+i = 0, and is easily proved by 
inserting fl4.8Q and using ( 4.10| ). 

As far as the last assertion is concerned, we use the following consequences of (|4.S|): 



(A* k p2k+l)+ = P2i+1, (A* k q2k+l)+ = <?2i+l • 



(4.11) 



This gives, using (fOl), 



^i-& v (2fc+l)) + _ V (2i+1) _ R 2i+l 





1 



Since —H^ t+1 is the coefficient of A fe 1 1 in t&ifc+ij equation ([4.6[) shows that we are done. 



yfc-i-1 : 



(4.12) 



□ 



So we have seen that the double reduction process of the Central System outlined in Section 
|H provides us with a natural Lax representation of the (commuting) vector fields of the KdV2 9 +i 
system. Actually, as it was explained in [15], the Central System can be seen as an outgrowth of 
the bi-Hamiltonian properties of the KdV hierarchy. It is thus natural to look for a bi-Hamiltonian 
structure of the KdV2 S +i system. Unfortunately, we are not in a position to derive such a property 
from the Central System itself, but rather we have to rely on the Lax representation discussed so far. 
Namely, in the next two sections we will establish the bi-Hamiltonian nature of KdV^p+i, showing 
that it comes from the general theory of bi-Hamiltonian systems defined on matrices depending 
polynomially on a parameter. 
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5 Lax Equations and bi-Hamiltonian Systems 



In the previous section we have associated with every point of the phase space M 2g +i of KdV2 9 +i 
a Lax matrix \A 29+1 ) , and we have seen that this matrix gives a Lax representation of the flows. 
To give these flows a bi-Hamiltonian formulation, we will address in this section a general problem, 
concerning the relation between Lax matrices and bi-Hamiltonian structures. We will describe a 
class of bi-Hamiltonian manifolds whose (bi-)Hamiltonian flows have a Lax formulation, and show 
that this formulation survives a reduction process of Marsden-Ratiu type. Since 

V(2fl+i) depends 

polynomially on A, it is quite natural to consider the multi-Hamiltonian structures defined on re- 



valued polynomials (see [33, 24]), where q is a Lie algebra of matrices such that the trace of the 
product is nondegenerate. More precisely, for a fixed matrix Agj, let us consider the space 



M A := {X{\) = X n+1 A + X * X i \XiEQ} 



(5.1) 



i=0 



which is clearly in a 1-1 correspondence with the space ©™ =o of (n + l)-tuples of matrices in q. 
The tangent and the cotangent space at a point of Ma can also be identified with ©™ =0 Q, using 
the pairing 



((V , . . . , V n ), (W , ...,W n )) = J2 Tr( W*) • 



(5.2) 



If F is a function on Ma, we will denote its differential by 

, , , , OF dF 
dr 



dXo'-'-'dXr, 

It is known that on Ma there is an (n + 2)-dimensional web of (compatible) Poisson brackets, and 
that this web is associated with a family of classical i?-matrices. Nevertheless, it turns out that in 
our case the relevant Poisson pair is given by the first two brackets of the above mentioned family. 
The first Poisson tensor, as a map from the cotangent to the tangent space, is given by 



Pn 



Wi 



\ 



I 



x 

Xi 



\ 



while the second one is 



X n -i 



X 
Xi 



( 



[Xi, 

[x 2 , 



y W n -! J y X n _! J \ [A, •] 



( - \XoA 








[x 3 , 



[A, 



[x 3 , ■] 



[A, 



\ 



( 



W 
Wt 



\ 



J \ Wn-i J 



\ / Wq ^ 



[A,-] 



Wi 



/ ^ Wn-i j 



(5.3) 



(5.4) 
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Figure 1: The bi-Hamiltonian reduction 



The associated Poisson brackets are given by {F, G}i = (dF,PidG), where i = 1,2, and F, G 
are functions on Ma- The Poisson tensors ( |5.3| ) and Q5.4| ) satisfy the remarkable property that 
every linear combination of them is still a Poisson tensor. For this reason one says that they are 
compatible, and that Ma is a bi-Hamiltonian manifold. 

Let us consider now the Poisson pencil P\ := P\ — XPq. It is important to notice that its 
Hamiltonian vector fields admit a Lax representation, as shown in the following: 

Proposition 5.1 Let F be a function on Ma, an d — » = P\dF the Hamiltonian vector field 

dt x 

associated by P\ to F. Then, 



dX(X) 
dt x 



dF 

-,X(X) 



dX 



(5.5) 



Proof. Use the expressions ( |5.3| ) and (5.4) to compute the vector field P\dF, then identify the 
parameter A appearing in the Poisson pencil with the A in X(X). 

□ 

At this point we want to enlarge the class of bi-Hamiltonian manifolds giving rise to systems 
admitting a Lax representation, having also in mind the case of the stationary reductions of KdV. To 
this aim, it is important to recall a reduction theorem || allowing us to "move" the bi-Hamiltonian 
structure from a given ( "big" ) manifold to a smaller one. This result is a particular case of a theorem 
by Marsden and Ratiu for Poisson manifolds [p6[ , and can be applied to a general bi-Hamiltonian 
manifold. Here, for the sake of simplicity, we will describe it only in the case at hand. The central 
point is that a Lax representation can be found also for the vector fields that are Hamiltonian with 
respect to the reduced Poisson pencil. 

The first step of the reduction process is to fix a symplectic leaf S of Po . Then we introduce the 
distribution D = Pi(KerPo)) which is integrable thanks to the compatibility between Po an d P\- 
From the explicit form ([5.3|)- (|5.4| ) of the Poisson tensors, it is easy to see that the vector fields in 
D have the Lax form X(X) = [Wo, X(X)], for a suitable Wq. Let us denote by E the intersection of 
D with TS. The statement of the bi-Hamiltonian reduction theorem is that the quotient manifold 
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TV" = S /E inherits from Ma a bi-Hamiltonian structure. In order to compute the reduced Poisson 
bracket {/, g}\ between two functions /, g on TV", we consider them as functions on S, invariant 
along the leaves of E. Then, we choose functions F and G on Ma which extend / and g, and 
annihilate the distribution D. Their Poisson bracket {F, G}\ is still invariant along D, and therefore 
defines a function on TV", which is independent of the choice of the prolongations F and G. 

Let us consider now a given Hamiltonian vector field Xf (with respect to the Poisson pencil) 
on TV", with Hamiltonian /. If F is a prolongation of /, the vector field Xp '■= P\dF is easily 
seen to be tangent to S and to project onto Xf. We are going to show that Xf inherits a Lax 
representation from the one of Xp. To this aim, we suppose that there exists a submanifold Q of 
S, which is transversal to the distribution E. In other words, Q is the image of a section of the 
bundle ir : S —* TV". Then Q is diffeomorphic to TV, and inherits a bi-Hamiltonian structure. The 
representative of Xf on Q is simply found by decomposing the restriction of Xp to 5 according to 
the splitting TS = TQ © E. Since we have seen that E is spanned by vector fields having a Lax 
form, we have proved the following: 

Proposition 5.2 Let Q C S be transversal to E. Then, the vector fields on TV" which are Hamil- 
tonian with respect to the Poisson pencil admit a Lax representation on Q. 

Therefore the bi-Hamiltonian reduction implies, in this case, a reduction of the Lax formulation. 



6 The bi-Hamiltonian Structure of the KdV Stationary Reduc- 
tions 

The aim of this section is to show that the KdV2 S +i systems introduced in Section || admit a 
bi-Hamiltonian formulation. To do this, we are going to exploit the Lax representation found in 
Section |4] and the results of the preceding section. 

The form of the Lax matrix \A 2ff+1 ) found in Section || suggests to choose q = s((2), n = g, and 



A 





1 



Therefore, the dimension of Ma is 3(g + 1). The Lax matrix \/^ 2 9 +1 ) defines an embedding of the 
KdV2 5 +i phase space into Ma- At this point two natural questions arise: 

1. Does this submanifold inherit from Ma the bi-Hamiltonian structure? 

2. If so, are the vector fields of KdV2 3 +i bi-Hamiltonian with respect to this structure? 

We will see that the answer to both questions is yes. 

In order to answer the first question we need a careful description of the symplectic leaves of 
Po, as given by 

Lemma 6.1 The symplectic leaves of Pq have dimension 2{g+\). Moreover, letH(X) : Ma — > C be 
defined as H(X(X)) := ^TrX(A) 2 and let the Hi be the coefficient of X 1 in H(X). Then, the functions 
H2 g +i, ■ ■ ■ , H g+ \ are functionally independent Casimirs of Pq. Consequently, the symplectic leaves 
of Pq are the level surfaces of the previous Casimirs. 

Proof. From ( |5.3| ) the kernel of Po is easily seen to be given by the covectors [Wo, . . . , W g ] T such 
that Wi = ctiA + Yl]=i a i-iXg+\-u where the a%, i = 0,...,g, are arbitrary. This shows that 
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dim(KerPo) = n, so that the dimension of the symplectic leaves is 2(g + 1). In order to check 
that H^g+i, ■ ■ ■ ,H g+ i are Casimirs, it is sufficient to verify that the differential of Hi is the 1-form 
[Xi, Xi-i, Xi-g] , where X k := if k < 0. 

□ 

One can easily show that the symplectic leaf defined by Hi = q, with g + l<i<2g + l, can 
be parametrized as 



X(X) = X 9+1 A + V 

3=0 



(6.1) 



where pj and qj are free parameters, and rj is a function of (j>j+i, qj+i, ■ ■ ■ ,Pg,Qg) and the values 
(cg+j+X, . . . , C2 5 +i) of the Casimirs. 

As far as the distribution D = Pi(KerPo) is concerned, in this case it is tangent to the symplectic 



leaves of Pq. Indeed, from the explicit form (|5.3;)-(5,4) of the Poisson tensors it is easy to see that 



D is the 1-dimensional distribution spanned by the vector field 

X(X) = [A,X(X)} . (6.2) 

This also shows that the integral leaves of D are simply the orbit of the action given by simultaneous 
conjugation of the isotropy subgroup of A, but we will never use this fact. 

Now we are ready to endow the phase space M.2g+i of KdV2 S +i with the structure of a bi- 
Hamiltonian manifold. This follows from the fact that the map assigning to each point of M.2 g +i 
the corresponding Lax matrix \A 2fl+1 ) defines a sub manifold of a suitable symplectic leaf of Pq, 
which is transversal to the distribution E. This is shown in the following: 

Proposition 6.2 Let us take the symplectic leaf S defined by #2g+i = 1 an d Hi = for g + 1 < 
i < 2g. Then, \/( 29+1 \hi, . . . , hag+i) G S for all (h\, . . . , /i2g+i) G A^g+l- Moreover, the image Q 
of the previous map is transversal to E. 

Proof. By the definition of S, we have to show that ^Tr (V^ +1 )) 2 = A 2fl+2 + Yl=o H i Xl - To 
this aim, we observe that equation (4.1) and the stationarity of ^29+1 imply that H^ 29+l \z) is 



an eigenvalue of \A 2ff+1 ). The other eigenvalue is given by #( 2 9 +1 )(-z), because V^ 1 ) depends 
only on A = z 2 . Therefore, Tr (V^ 1 )) 2 = (H^ +1 )(z)) 2 + {H^ +1 ) (- z)) 2 has the desired form. 
Finally, referring to the parametrization ( |6.lD , one can easily prove that the submanifold p g = is 
transversal to the distribution E. 

□ 

Hence, the KdV2 9 +i phase space inherits from Q (and from the quotient space M = S/E) a bi- 
Hamiltonian structure. To compute this structure, it is convenient to use the formalism discussed 



in [10], whose aim is to avoid dealing with the explicit form of the projection ir : S — > M. 

Now we want to show that the KdV2 9 +i flows are indeed bi-Hamiltonian with respect to the 
above Poisson pencil. The Lax representation ( |4.4j ) and the form ( |4.5| ) of the Lax pair suggest to 
consider the vector fields on Ma given by 

dX(X) 



X l ~°X(\)) X(X) , i = -1,0,... ,</-!. (6.3) 



They are Hamiltonian with respect to the Poisson pencil on Ma, with Hamiltonian function given 
by X(X) ^ (X l -3H(X{X))) where H{X(X)) = ^TrX(X) 2 . Furthermore, we can state 
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Proposition 6.3 On the bi-Hamiltonian manifold Ma the function H(X) ■ 
Casimir of the Poisson pencil. The bi-Hamiltonian vector field Y g —i ■= PodHi 



eEo 1 m « 

= P\dHi has the 



Lax representation (6.8 ). 



Proof. One can easily see that dH{\) = [X(X), XX (X), . . . , X 9 X(X)Y . Thus from (fj) and 
Q5.4j ) it follows that P\dH(X) = 0. The vector fields Yi are Hamiltonian with respect to P\, since 
Yg-i = P\dHi = P\d(X~ l H(X)) + . Thus the Lax representation (6.3) is a consequence of Proposition 

□ 



Remark 6.4 The previous proposition is a particular case of a general result |33|], stating that 
Ad-invariant polynomial functions on a Lie algebra q give rise to Casimir of the Poisson pencil on 
M A . 



In order to obtain the KdV2 9 +i system from the bi-Hamiltonian vector fields (p. 3]), we remark that, 
from general results of the bi-Hamiltonian theory: 

1. The functions Hi are invariant along the distribution E, and therefore can be projected on 
the quotient J\f. 



2. The vector fields ( |6.3[) are tangent to S and project onto M. 

3. Their projections are the bi-Hamiltonian vector fields associated with the projected functions. 



We observe from ( |6. 3[ ) that the vector field Y—i = PodH g is tangent to the distribution E. This 
means that the function H g , on the quotient TV, is a Casimir of the reduction of Pq. Hence, the 
polynomial Hq + H\X + • • • + H g X 9 is a Casimir of the reduced Poisson pencil. The other vector 



fields in (6.3) project on the stationary reductions of KdV, as shown in 



Proposition 6.5 The projections on J\f of the vector fields (6.5), for i = 0, . . . ,g — 1, coincide 
with the KdV2 g +i systems. 

Proof. The right place to compare the two hierarchies of vector fields is the transversal submanifold 
Q defined as the image of the map V( 2fil+1 ) : M2 g +i — ► <S. Hence, we must project the Lax equations 
( |6.3| ) on TQ along E. This leads to 

(A^X(A)) X(A)1 - a [A,X(X)} , 



(6.4) 

his means that the en- 
{X^X(X)) +1 X(X)~ = 

" 4_1 -coefficient of X(X), 
□ 

Therefore, we have shown that the KdV2 9 +i flows are bi-Hamiltonian. Moreover, they are 
associated with a Casimir of the Poisson pencil, having a polynomial dependence on A. Since this 



where a is fixed by the condition that this vector be tangent to Q. 
try (1, 1) of the coefficient of A 9 in ( |6.4[ ) must be zero. If we write 

X(A),(A^X(A)). 



, then we obtain that a is the entry (1,2) of the X 9 ' 



so that equation (4.5) concludes the proof. 



is a particular instance of a general theory developed by Gel'fand and Zakharevich [19|, we will say 
that the vector fields of KdV2 9 +i are GZ systems. The next section is devoted to the separability 
of such systems. 
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7 Separability of bi-Hamiltonian Systems 



We have just seen that the stationary reductions of KdV are examples of GZ systems. In this section 
we show that the bi-Hamiltonian structure of such systems allows one to solve them by separation 
of variables. Under special circumstances, separability of GZ systems was proven in |5|, |2?J. We 



refer to jl6| for complete proofs and a more detailed discussion. 

Let A4 be a (2n + l)-dimensional manifold endowed with a pencil P\ = P\ — XPq of Poisson 
tensors. We suppose that the rank of P\ is generically 2n, so that (locally) there exists a polynomial 
Casimir function H(X) = Y17=o °^ ^ ( see EH)- The associated GZ systems PodHi = P\dHi + i 
are obviously tangent to the symplectic leaves of Pq, and give rise to Liouville integrable systems. 
Since H n is a Casimir of Pq, such leaves are the level surface of H n . Let us denote by u the 
symplectic form given by the restriction of Po to a (fixed) symplectic leaf; if Xf := Podf, where / 
is any function on M, then 

u(X f ,X g ) = {f,g} . 

In order to exploit the existence of the other Poisson bracket, we make an additional assumption. 
We suppose that there exists a vector field Z on M. such that 

1. It is transversal to the symplectic leaves of Pq; 

2. The functions invariant along Z form a Poisson subalgebra with respect to the bracket {•, -}\ 
associated with the Poisson pencil. 

The second condition means that the bi-Hamiltonian structure can be projected on the quotient 
space of the integral leaves of Z. The first one tells us that such quotient can be identified with a 
symplectic leaf S c of Pq, which is therefore a bi-Hamiltonian manifold. Moreover, we can define on 

S c a Nijenhuis tensor N as 

u(X f ,NX g ) = {f,g} 1 , 

where / and g are functions on Ai invariant along Z. Thus S c is said to be a Poisson-Nijenhuis 



(PN) manifold (see [21] and references cited therein). 

The vector field Z allows us to use the Poisson pencil to construct variables of separation 
for the GZ systems on <S C . Hence in this case the Poisson pencil not only provides us with a 
commuting family of vector fields, but also gives coordinates for which the corresponding equations 
of motion can be solved by separation of variables. Indeed, one can show that the Nijenhuis 
tensor N has n functionally independent eigenvalues (Ai,...,A n ). Then (see, e.g., f25f|), there 
exist n complementary coordinates (fj,\, . . . , fj, n ) on S c such that to takes the canonical form to = 
Y^i=x d\ A d/ii and the adjoint N* of N takes the diagonal form 

N*dXj = AjdAj, N*dfij = Xjd^j . 

Such coordinates are called Darboux-Nijenhuis (DN) coordinates. They are the separating coor- 
dinates for the GZ systems. Indeed, let us normalize Z in such a way that Z{H n ) = 1. Then 
the differentials of the restrictions Hi of the Hamiltonians Hi to S c generate a subspace which is 
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invariant with respect to N*. More precisely, we have that 



N*dH, 



o 



••• c 

1 ••• ci 

1- C 2 







1 c n _i 



n-l 



(7.1) 



where Cj = —Z(Hi). This implies also that 

minimal polynomial of N = A n - ^ qA* = Z(H(\)) . 



n-l 



(7.2) 



j=0 



Moreover, one can check that the Frobenius matrix F defined by ([O]) satisfies the condition 

N*dV = F <iF , (7.3) 

where dF is the matrix whose entries are the differentials of the entries of F, and, on the left-hand 
side, N* acts separately on each entry. Conditions (|7.3|) and (7.1) imply that the Hamilton- Jacobi 



equations for the Hi are (collectively) separable in the DN coordinates. In fact, the (transpose of 



the) Vandermonde matrix constructed with the \j diagonalizes F, and applied to [Hq, 
gives a Stdckel vector, in the sense that its j-th component depends only on (Xj,/jLj). 

As far as the explicit construction of the DN coordinates is concerned, we have seen that the Xj 
are the roots of the minimal polynomial Z{H{\)) of N. On the contrary, the coordinates /i,- must 
be computed (in general) by a method involving quadratures. However, in the case at hand there 
is a recipe that is particularly useful in the applications. Let us consider the Hamiltonian vector 
field Y on S c , associated with ^TVA r = X^i^i by the symplectic form u>. If the (restriction of the) 
Casimir H(X) satisfies the condition Y r (H(X)) = for some r, then the coordinates 



f'.i 



Y r ~ 1 (H(Xj)) 



form with the Aj a set of DN coordinates. Hence in this case the bi-Hamiltonian structure provides 
us with a method to algebraically construct the separation variables. 



8 Separability of the Stationary Reductions 

In this section we will show that the KdV2 9 +i system belongs to the class of separable GZ systems 
discussed above. It is convenient to show that the conditions on the vector fields Z and Y are 
fulfilled on the "big" bi-Hamiltonian manifold Ma and then to reduce everything. 

Regarding the transversal vector field Z, we introduce on Ma the vector field Z Ma defined as 

X = A, Xi = for all * = 1, . . . , g. (8.1) 

It is tangent to the symplectic leaves of the Poisson tensor P$ given by (|5.3|) (since it is easily seen 
to belong to its image) and it can be projected on the quotient space TV = S / E (since it commutes 
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with the generator (|6.2| ) of the distribution E). Using again the form (|5.3| ) and ( |5.4| ) of the Poisson 
tensors on A4a, one can check that the functions invariant along Z Ma form a Poisson subalgebra 
with respect to P\. This property is trivially conserved after the reduction on J\f. Finally, we have 
to show that the reduced vector field is transversal to the symplectic leaves of (the reduction of) 
Pq. This follows from the fact that, at the points where Tr(X 9 ^4) = 1, 

L z m a (iTrX(A) 2 ) = Tr (X(X)A) = X 9 -\ , 

so that L Z M A H g = 1. This also shows that the reduction of Z Ma has the right normalization. 

Thus, we have shown that on Af, which is diffeomorphic to the phase space of KdV2 9 +i, there 
exists a vector field that satisfies the hypotheses of the previous section. Therefore, the stationary 
reductions of KdV can be solved by separation of variables in the DN coordinates. We are left with 
the problem of finding explicitly these coordinates. 

To this aim, we introduce the vector field Y Ma defined on Ma &s 

X = [A, X g ], X i = for alH = 1, . . . , g. (8.2) 

It is also tangent to the symplectic leaves of the Poisson tensor Pq and can be projected on the 
quotient space N . Moreover, 

1. Y Ma is (up to a sign) the Hamiltonian vector field associated by means of Pq with the Lie 
derivative along Z Ma of the coefficient Hg-i of ^TrX(A) 2 ; 

2. We have that L\ Ma (H{X)) = 2 (Tr(AX g )) 2 . 

The first assertion can be checked after noticing that Z MA (H g -.\) = Tt(X g -iA), and that the 
differential of this function is [0, . . . , 0, A, 0] T . The second assertion simply follows from the fact 
that H(X) = ±TrX(A) 2 . 

The same properties hold also on Kf: If Z and Y are the reductions of Z Ma and Y- Ma , we 
have that Y is the Hamiltonian vector field associated with —Z(H g -i) by the reduction of Pq. 
Furthermore, we have that Ly{H{X)) = 2, since Tr(AX g ) = 1 at the points of S. Let us now 
restrict ourselves to a symplectic leaf of the reduction of Pq. Since —Z(H g _i) = c g -\ = ^TriV, the 
vector field Y can be used to construct the pLj coordinates according to the recipe given at the end 
of the previous section. The conclusion is: 

1. The Xj are the roots of the polynomial Z{H(X)) = Ti(AX(X)) = Tr^V^ 1 )), that is, the 
entry (1,2) of the matrix \A 29+1 ); 

2. The fij are given by fij = f(Xj), where 

/(A) = y^hm) = l T "( x ( x ^ x ^ = entr y ( 2 > 2 ) of v(29+1) - 



We remark that our general theory of separability gives, in this particular case, the same construc- 
tion of the variables of separation holding for systems admitting Lax with parameter formulation 
(see, e.g., Q |3(|). In fact, writing the Lax matrix (|4.6| ) as 



VC29+ 1 ) 



Vg(X) Ug{X) 
Wg{X) -Vg{X) 
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the equation of the associated spectral curve C is 

f = U g (X)W g (X) + V g (X) 2 . 

Since U g (\j) = and fij = —V g (Xj), we see that g points (Xj,fXj) lie on C. 

We close this section with a description, from our point of view, of the example of KdVs we 
started with in Section 0. We consider the bi-Hamiltonian manifold 



M A 



AX 3 + X 2 X 2 + XxX + X \X t 



Pi 



-Pi 



whose Poisson tensors are given by (5J3) and (5.4). The reduction process described in Section || 
allows us to pass to the transversal submanifold 



,4A 3 + 


1 


A 2 + 


Pi 


-<?2 


x + 


Po Q2 2 - qi 




Q2 






-P\ 




qo -po 



which is diffeomorphic to the phase space M§ of KdVs. The correspondence is given through the 



Lax matrix \A 5 ) displayed in 


Section ^ 


J. The resulting chang 


e of variables is explicitly 




h = q 2 , 


h 2 = 


-Pi , h 3 = qi , 


hi = -p - piq 2 , 








h 5 = 


<?i<?2 + \pi - l-q2 3 + 


ho ■ 


The Poisson pencil 


on M.5 turns out to be 









-1 





—hi + A 


-h 2 




1 





2hi - X 


h 2 


h 3 + \hi 2 - 2hiX 


Px = 





-2hi + X 





-ha - hi 2 + 2hiX 


—hi — hih 2 




hi - X 


-h 2 


* 





-h 5 + 3hih 3 - \h 2 2 - hi 3 








(2h 3 + hi 2 )X 




* 


* 


* 


* 






(8.3) 



Its Casimir H(X) 
matrix: 



Hq + HiX + H 2 X 2 can be computed with the trace of the square of the Lax 



2h 3 h 5 + hi 5 + 2hih 3 2 - 2hih 2 h A - Zh 3 h 3 + 2h 2 h^ + h 4 2 



Ho = h 3 h 2 
Hi = 2h 2 h 4 - 2hih 5 + 3hi 2 h 3 
H 2 = 2hi 3 - Ah x h 3 + 2h 5 



hih 2 2 



h 3 2 



hi* 



The two vector fields of the KdVs hierarchy are given by (3.12). The symplectic leaves of Po are the 
level surfaces of H 2 . The vector field Z Ma is d/dqo, while its projection Z on M.§ is (\/2)d/dh§. 
On the symplectic leaf S c defined by H 2 = c we can use (hi, h 2 ,h 3 , h±) as global coordinates, and 
the corresponding Poisson pencil is simply obtained by deleting the last row and the last column 
in (8.3). The minimal polynomial of the Nijenhuis tensor on S c is 



Z(H(X)) = X 2 - hiX + h\ - h 3 
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and Ai, A2 are its roots. To find fi± and [i 2 we have to use 

yM-A . 



d d 

r2 a h 2 P 2 a~ ' 

opo oq 



whose reduction on M 5 is Y = d/dh 4 . Since Y(H(X)) = 2h 2 X + 2h 4 - 2h 1 h 2 and Y 2 (H(X)) = 2, 
the coordinates fj,\ and [i 2 are the values of the polynomial 

h 2 \ + /14 — h\h 2 

for A = Ai, X 2 . In order to check that the DN coordinates are separation variables for the restrictions 
Hq and Hi of the Hamiltonians to S c , we simply have to compute 



1 Ai 




H 




-cAi 2 + ijli 2 - Ai 5 


1 A 2 




. Hi _ 




-cA 2 2 + ijl 2 2 - A 2 5 



from which the form of the spectral curve can also be seen. 



9 Final Remarks 



1. The results outlined in Section are proved in [16] for a class of bi-Hamiltonian manifolds whose 
rank is not maximal. This means that our approach to the stationary reductions of KdV can be 
directly generalized to the stationary reductions of the Gel'fand-Dickey hierarchies. A step in this 



direction has already been taken in [17], whose results should be compared with those of IS, 31]. 
We will treat this problem in a future publication. 

2. The separation variables provided by the bi-Hamiltonian method coincide, in the KdV case, 
with the ones obtained by algebro-geometric constructions. It would be interesting to compare in 



more general cases these two methods. A first result has been obtained in 2C], where the "spectral 
Darboux coordinates" of ^ are shown to be DN coordinates for a suitable pair of compatible 
Poisson brackets. 

3. Another Marsden-Ratiu reduction of the manifold Ma of Section || has been performed in [30], 
for an arbitrary simple Lie algebra g. That reduction leads to a bigger quotient space, and allows 
one to reduce all the multi-Hamiltonian structure of Ma and to obtain, in the case g = sl(2), the 
Mumford systems [25]. A further restriction to the level surface of some Casimirs gives the same 
reduced phase space obtained in Section |(| where only two Poisson brackets survive. 
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